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1.  INTRODUCTION 


f 


% 


This  report  is  concerned  with  results  of  the  first  part  of  a feasibility 
study  of  an  improved  air  filtration  system.  It  uses  ionization-producing  beta- 
ray  isotope  radiation,  in  conjunction  with  electric  charging  and  acceleration, 
to  cauge  fine  particles  to  agglomerate,  hence  enhancing  the  subsequent  collec- 
tion process. 

The  removal  of  aerosol  particles  from  a gas  by  any  of  the  standard  gas- 
cleaning techniques  (filtration,  scrubbing,  or  electrostatic  precipitation)  is  known 
to  be  more  difficult  for  submicron  particles  than  for  larger  particles.  Filtration 
methods  presently  in  use  can  achieve  high  efficiencies  of  submicron  particle  re- 
moval, using  high- efficiency  HEPA  filters.  These  filters  are  costly  to  build 
and  replace.  Moreover,  frequent  replacement  as  is  necessary  in  some  applica- 
tions implies  a large  amount  of  system  shutdown  time.  In  the  present  state-of- 
the-art  inexpensive  low- efficiency  prefilters  are  used  ahead  of  the  HEPA  filters 
to  capture  the  bulk  of  the  aerosol  mass  residing  in  the  larger  particles.  The 
finer  particles  not  captured  by  the  prefilters  then  pass  into  the  HEPA  filters. 

To  the  extent  that  the  prefilters  can  be  made  more  efficient  in  capturing  fine 
particles,  the  lifetime  of  the  HEPA  filters  can  be  prolonged,  with  a correspond- 
ing reduction  in  cost  and  shutdown  time. 

The  system  investigated  is  intended  to  agglomerate  fine  particles  into 
larger  sizes  so  that  they  are  more  easily  captured  by  the  prefilters. 

The  following  sections  give  details  of  the  filtration  enhancement  system 
andthe  key  physical  processes  including  ionization,  particle-charging,  andagglom- 
eration.  The  physical  parameters  of  the  system  are  presented  first;  then  the  elec' 
trical  and  operational  parameters  are  determined  which  are  required  for  the  sys- 
tem to  demonstrate  significant  radiation-induced  particulate  agglomeration.  A 
preliminary  evaluation  of  the  system  feasibility  is  then  given;  the  report  concludes 
with  recommendations  relative  to  the  further  experimental  and  analytical  work 
necessary  to  fabricate  a model  and  to  verify  the  predicted  filtration- enhancement 
characteristics  of  the  system. 


2.  the  filtration  enhancement  system 

2.1  Basic  Description 

The  proposed  system  consists  of  a charging  stage,  followed  by 

one  or  more  agglomerating  stages  (see  Fig. 2.1).  As  illustrated,  each 

stage  is  essentially  a cylindrical  tube.  Within  each  tube  the  gas  flows 

along  isolated  sections,  of  pie-shaped  cross  section,  partitioned  from 

one  another  by  radial  vanes  (see  Fig. 2. 2).  Along  the  axis  of  the  first 

tube,  which  is  the  charging  stage,  is  an  inner  cylinder  containing  radio- 
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active  material  (for  example,  Kr  gas).  In  this  stage,  the  gas  outside 
the  inner  tube  is  ionized  by  the  radiation  passing  through  the  wall  of  the 
inner  tube.  DC  voltages  of  alternating  sign  are  applied  to  alternating 
partition  vanes  (Fig. 2.2).  These  applied  voltages  cause  the  ions  to  con- 
centrate in  separate  unipolar  regions.  The  gas  containing  the  particles 
to  be  agglomerated  flows  through  these  unipolar  regions  by  a flow- 
splitting geometry  (Fig.  2.2). 

After  the  particles  become  charged  in  the  charging  stage,  they 
pass  into  the  first  agglomerating  stage.  This  stage  has  partitioning 
vanes  identical  to  and  aligned  with  those  of  the  first  stag-?,  but  the  signs 
of  the  applied  voltages  are  reversed.  This  forces  oppositely -charged 
particles  to  move  toward  each  other  so  that  they  would  tend  to  agglomerate 
by  collision.  Particles  failing  to  agglomerate  in  the  first  agglomerating 
stage  pass  into  the  next  agglomerating  stage  where  the  vane  voltages  are 
again  reversed.  Hence,  the  particles,  once  charged,  are  accelerated 
toward  oppositely-charged  partners  in  successive  passes  in  a sequence 
of  agglomerating  stages.  The  agglomerated  particles  can  be  separated 
from  the  unagglomerated  particles  by  an  appropriately  designed  (e.  g.  , 
impactor  type)  gas -flow  system. 


Ionization- Induced  Charger 
showing  ion  distribution  in 
typical  section 


Typical  diameter  50-100  cm 


Radioactive  gas  in  thin-walled 
container  to  allow  betas  to 
escape  with  minimum  energy 
loss 


Vanes  with  DC  voltage 
of  alternating  sign 


Collimator  to  limit 
region  of  ionization 
to  center  of  pie- shaped 
s ection 

Ions  are  drawn  toward 
plate  of  opposite  sign 


Beta  particles  deposit  / 
energy  in  ionization  region 
producing  charged  ions 


Agglomeration 

region 


After  leaving  Charger, 
charged  particulates  * 
enter  near  vane  of  same 
sign  and  move  toward 
particles  of  opposite 
charge 


* ""  Reverted- Field 

Agglomerator 

Fig.  2.2.  Cross  Sections  of  Charger  and  Agglomerator  Showing  Typical  Pie-Shaped  Section 


2.2  Areas  of  Investigation 

As  discussed  above,  the  basic  concept  involves  use  of  a radiation 
source  to  produce  separate  regions  of  positive  and  negative  ions  through 
which  the  particulate- containing  air  passes  (Fig.  2.2).  After  leaving  the 
charging  stage  the  oppositely-charged  particulates  enter  an  agglomeration 
stage  (Fig.  2.1)  in  which  the  particles  of  opposite  sign  are  brought  into  close 
proximity,  hence  producing  agglomeration  and  shifting  the  particle  distribu- 
tion toward  larger  sizes.  The  result  is  that  the  average  efficiency  of  the 
following  collection  system  - be  it  scrubber,  prefilter,  etc.  - is  increased 
provided  the  average  particle  size  is  shifted  past  the  well-known  minimum 
(near  0.2-0.  3f/m  particle  size)  in  the  final  stage  collection  system.  The 
geometry  for  implementation  of  this  concept  as  shown  in  Figs.  2.1  and  2.2 
should  be  considered  as  illustrative  only,  not  limiting.  Other  geometries 
(for  example,  rectangular)  involving  the  same  concept  can  be  imagined. 

There  are  five  general  areas  requiring  investigation: 

1)  Production  of  ionization  by  the  radiation  source. 

2)  Separation  of  ionization  to  produce  a distribution  having 
regions  of  predominantly  positive  and  negative  ions. 

3)  Charging  of  the  particulates  by  the  ions. 

4)  Agglomeration  of  the  charged  particles  to  shift  distribution 
toward  larger  sizes. 

5)  Collection  (filtration)  of  the  altered  particulate  distribution. 

As  a first  step  in  demonstrating  feasibility  of  this  concept  it  is 

necessary  to  show  that  sufficient  ionization  can  be  generated  by  realistic 
activity  radioactive  sources.  The  possibility  of  using  ^Co  gamma 
radiation  to  generate  ionization-producing  electrons  in  a duct  wall  has 
been  investigated  by  Heinsohn  et  al.  (Ref. 2.1).  Their  intent  was  to  pro- 
duce charged  particulates  for  purposes  of  electrostatic  collection,  a 
completely  different  objective  fro?r  the  present  work.  Furthermore, 
they  did  not  separate  the  region  of  ionization  production  from  those  in 
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which  the  positive  and  negative  ion  distributions  are  maintained,  as  is  done 
in  Fig.  2.2.  We  believe  that  a beta  source,  which  emits  electrons  directly, 
will  prove  to  be  of  greater  utility  from  the  point-of-view  of  both  handling 
and  ionization  production,  than  will  a high  energy  gamma  source.  Krypton-85 
appears  to  us  to  be  of  particular  interest  because  of  its  availability,  specific 
activity,  and  the  range  of  its  betas  (up  to  ~200  cm).  Hence  the  basic  feasi- 
bility  of  producing  ionization  with  this  source  is  considered  in  some  detail 
below.  This  is  followed  by  a more  general  discussion  of  the  other  areas 
above,  including  the  analytical  methods  applied  to  each. 

In  order  to  carry  out  the  investigations  of  charging  and  agglomeration, 
it  is  desirable  to  consider  particular  particle  size  distributions.  A number 
of  typical  particle  size  distributions  encountered  in  nuclear  materials  hand- 
ling facilities,  and  characterized  as  log-normal,  are  summarized  in  Table 
2.1.  This  table  is  taken  directly  from  a report  by  Lawrence  Livermore 
Laboratory,  where  enhanced  filtration  mechanisms  of  a different  type  than 
that  of  interest  here,  are  being  investigated.  It  is  seen  that  the  particulates 
generally  cover  the  range  of  about  0.2  to  30  microns,  with  median  diameters 
of  a few  microns  and  standard  deviations  of  2-5. 
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Table  2.1 


Typical  particle  distributions  encountered  in  nuclear  materials 
handling  facilities.  (From  Ref.  2.2) 


0 

Fac  1 1 1 ty 

Process 

Aerodynamic 
mass  median 
diameter  (i.nl 

Ceomotr Ic 
std.  dev. 

(V 

Range 

(bm) 

Principal 

constituents 

R&D 

1.8 

2 . 2 

0.4  - 8.7 

Pu 

Ril) 

2.1 

3.3 

0.2  - 28 

Pu 

Kabr lc.it  ii 

>n 

4.0 

1 . 7 

1.4  - 12 

Pu 

Recovery 

0.  34 

3.4 

0.01  - 9.9 

Pu 

Kabrlc.it  ion 

2.  7 

2.4 

0.5  - 16 

Pu 

Hooky  Flats^^ 

Recovery 

0.9a 

vj 

0 

u> 

1 

Pu,  HF,  HNO j 

"Demonst rat  Ion 

. ' 

Planl"b 

Grinding, 

730"C 

11.9 

1.39 

- 

Pu  U 

"Industrial  Plant"*1 

Grinding. 

730°i: 

1.88 

1.76 

- 

Pu  U 

"Demons t rat  ion . 

Industrial  Plant"*3 

Grinding, 

1 750°C 

2.  30 

1.65 

- 

Pu  l! 

"Industrial  Plant"*1 

Grinding, 

I750"r. 

2.09 

1 .61 

- 

Pu  U 

3 Count  me d lan  diameter. 

**  G.  Newton,  Lovelace  Foundation,  private  conunictt Ion. 


References  6)  and  7)  above  are  as  follows: 

6)  H.  J.  Ettinger,  J.  C.  Elder  and  M,  Gonzales , Performance  of  Multiple 
HEPA  Filtere  Against  Plutonium  Aerosols,  LASL,  Report  No.  LA-5349- 
PR(1 973). 

7)  R.  W.  Woodward,  K.  J.  Grossaint,  and  J.  L.  McFeeters,  Exhaust  Fil- 
tration on  Gloveboxes  used  for  Aqueous  Processing  of  Plutonium,  14th 
ERDA  Air  Cleaning  Conference,  Sun  Valley,  Idaho  (1976). 
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3.  CYLINDRICAL  BETA  SOURCE  INDUCED  IONIZATION  DISTRIBUTION 


* * 


In  this  section  the  ionization  distribution  results  obtained  from  the 
85 

Kr  beta  source  are  presented  and  discussed  first.  This  is  followed  ty 
an  analytical  estimate  of  the  ionization  based  on  neglect  of  all  scattering 
effects.  This  is  known  to  be  a poor  approximation  if  the  actual  number  dis- 
tribution of  the  electrons  is  of  interest.  However,  it  is  found  here  that  rea- 
sonably  good  agreement  with  the  experimental  measurements  can  be  obtained 
by  this  simple  approach  when  only  the  ionization  distribution  is  of  interest. 

3.1  Experimental  Results 

85 

The  basic  configuration  of  the  Kr  gaseous  source  used  in  this 
work  is  shown  in  the  top  portion  of  Fig.  2.2.  The  cylindrical  beta  source 
itself  is  surrounded  by  a collimator  designed  to  limit  the  ionization  region 
in  the  manner  shown  there.  Fig.  3.1  is  a reduction  of  a full-scale  drawing 
of  the  source.  The  fabrication  is  completed  by  soldering  a thin  brass  foil 
(~2  mil)  over  the  entire  outer  surface  to  seal  the  assembly.  The  slots  are 
the  regions  through  which  the  betas  enter  the  collimator  apertures  in  Fig.2.2; 
the  solid  portions  of  the  source  would  be  obscured  by  the  collimator  walls. 
This  is  necessary  due  to  the  scattering  of  the  betas  after  passing  through 
foil. 


3. 1 . 1 Measurement  of  Source  Activity 

The  source  was  filled  with  gas  to  0.10  atm.  Since  the  specific  activity 
of  the  gas  was  not  precisely  known,  the  resulting  activity  of  the  source  was 
measured  with  the  method  outlined  in  Fig.  3.2. 


r 


/ 


8 


1/8  inch  A1  beta  shield 


Source  Activity. 


A well- shielded  (2  inches  of  lead  in  all  directions)  3 inch  diameterx  3 inch 
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deep  Nal(Tl)  crystal  was  used  to  detect  the  0.514  MeV  y-ray  emitted  by  Kr  in 
0.41%  of  all  decays.  A central  1/2  inch  diameter  portion  of  the  Nal(Tl)  crystal 
was  irradiated  by  a 3/4  inch  section  of  the  Kr  source,  all  other  portions  of  the 
source  being  shielded  by  2 inches  of  lead.  A piece  of  1/8  inch  aluminum  was 
used  to  absorb  the  emitted  beta  particles  and  shield  the  Nal(Tl)  crystal  from 
this  interference. 

The  total  source  activity  can  be  written  as 


S 

a 


€ F F F t. , i 
p o y s A1  s 


(3.1) 


where  S is  the  source  activity  in  disintegrations/sec,  C is  the  measured 
a p 

Nal(Tl)  count  rate  in  the  full  energy  peak  (0.514  MeV)  in  cps,  €p»0.80  is  the 
detection  efficiency  for  the  full  energy  peak,  F^  = 1.09x10"^  is  the  fractional 
solid  angle  (of  4^)  of  the  1/2  inch  diameter  collimator,  F = 4.1x10"^  is  the 


I 


o 
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yield  of  0.514  MeV  y-rays  from  Kr,  F = 6.4x10  is  the  fraction  of  the  source 

s 


viewed  by  the  Nal(Tl)  crystal,  tA1=  0.93  is  the  transmission  of  1/8  inch  A1  for  the 


0.514  MeV  y-rays,  and  t = 0.91  is  the  transmission  of  the  source  walls  for  the 


0 514  MeV  y- ray.  The  measurements  gave  S = 1.  04x1 01  °/sec,  or  S =0.28  Ci  of 
85Kt. 
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The  method  was  checked  by  also  measuring  the  activity  of  a 10//Ci  1 Cs 


calibration  source,  detecting  the  0.662  MeV  y- ray  emitted  in  85%  of  all  decays 


(Fy  - 0.85).  HereFg=l,  tg=l,  tA1=0.86,  f^arO.75,  and  Fq  is  the  same.  The  re- 


sult gave  10.4j*Ci,  and  thus  demonstrated  the  accuracy  of  the  measurement. 

85. 


The  measurements  show  that  the  gas  in  the  source  is  1.7%  Kr  by  volume. 
Since  the  active  volume  of  the  source  is  30.0  cm  long,  the  specific  activity  of 
the  source  is  0.0093  Ci/cm. 


3.1.2  Measurement  of  Ionization  Distribution 


The  Kr  ionization  source  was  placed  in  the  center  of  a four  ft.  long  by 
17-1/4  inch  ID  steel  pipe  (3/8  inch  wall)  and  the  radial  ionization  profile  mea- 
sured. This  geometry  was  designed  to  simulate  that  of  Fig.  2.2,  except  that 
the  vanes  and  collimator  are  not  present.  The  method  is  illustrated  below. 


Steel  pipe  4 ft.  long 
17-1/4"  ID 
3/8"  wall 


85 


Kr  source 
(center  ed) 


Sensitive  volume  of 
open-gridded  air 
ionization  chamber 
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Fig.  3.3  Setup  for  Measuring  Kr  Source-induced  Ionization 
Distribution. 
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The  ionization  chamber  is  of  the  open-gridded  air  type,  as  described 
in  Ref.  3.1.  An  outer  shield  mesh  is  held  at  ground  potential,  with  an  inner 
bias  mesh  1.5  cm  in  diameter  and  12.7  cm  long  defining  the  sensitive  volume. 
The  central  wire  is  grounded  through  an  electrometer  to  measure  the  ion  cur- 
rent, while  the  inner  mesh  is  held  at  + 100V  through  a 10  resistor.  Tests 
of  the  ion  chamber  showed  negligible  current  change  for  biases  more  than 
+ 30 V,  and  the  sign  of  the  bias  voltage  only  affected  the  sign  of  the  current 
measured,  not  its  magnitude.  The  shield  mesh  voltage  had  a negligible  ef- 
fect on  the  measured  current  in  the  absence  of  significant  air  currents,  and 

so  was  kept  grounded.  In  the  absence  of  ionizing  radiation  the  leakage  current 

-14 

was  measured  to  be  less  than  +1x10  A for  bias  grid  voltages  of  +100V.  The 
sensitive  volume  of  the  ion  chamber  is  22.0  cm3,  and  it  has  been  corrected 
for  the  volume  of  the  central  wire  and  bias  grid  support  structure. 

The  ion  chamber  was  scanned  over  center-to-center  distances  (radial, 

Kr  source  to  ion  chamber)  of  R=  5.5  to  17  cm,  and  the  current  measured.  The 
current  ranged  from  4.50xlO‘10A  to  1.26xlO_10A.  The  ion  production  rates  are 

Q = 1/  (eV.)  (3.2) 

where  Q is  in  ion  pairs/(cm  -sec),  I is  the  measured  ion  chamber  current  in 

- 1 9 3 

A,  e=  1.60x10  C is  the  electron  charge,  and  V.=  22. 0 cm  is  the  ion  cham- 

i 

ber  sensitive  volume.  The  resulting  Q profile,  normalized  to  the  source  inten- 
sity, is  shown  in  Fig.  3.4.  The  estimated  uncertainty  is  jt20%. 

At  distances  of  6-18  cm  from  the  source,  measurements  were  also  made 
of  the  ionization  distribution  along  a circumference.  Very  little  such  depend- 
ence was  found,  indicating  that  atmospheric  scattering  is  a large  effect,  par- 
ticularly in  the  absence  of  the  collimators  of  Fig.  2.2. 
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Specific  Ioniation  Rate  Q/S„  flO8  ion 
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Analytical  Result 
Neglecting  Scattering 


Experimental 

Data 


Outer  Edge 
of  Source 


Edge  of  Pip 


Source- ion  chamber  distance  R(cm) 

Fig.  3.4  Measured  and  Calculated  Specific  Ionization 
Distribution  from  8*Kr  Source. 
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3.2  Analytical  Determination 

Exact  calculation  of  ionization  distributions  around  sources  such  as 
that  of  Fig.  3.1  can  only  be  made  by  Monte  Carlo  methods,  in  which  scattering 
is  included.  Here  we  will  make  a calculation  based  on  the  approximation  in 
which  scattering  is  entirely  neglected. 

Consider  the  geometry  shown  in  Fig.  3.5,  which  shows  the  source  cross 
section  and  two  points  (a  distance  R from  the  source  center)  at  which  the  calcu- 
lation of  ionization  is  to  be  made.  Neglecting  scattering,  in  each  case  only  beta 
particles  emitted  from  gas  in  the  shaded  areas  can  reach  the  points  of  interest. 
Note  that,  for  the  distance  R shown  (equivalent  to  10  cm  in  actual  scale)  the  ef- 
fective area  is  about  the  same  in  each  case.  The  measurements  reported  in  3.1 
above  found,  in  fact,  that  there  was  little  dependence  of  the  ionization  rate  on 
circumferential  location.  Hence,  to  do  the  calculation  here  we  can  replace  the 
actual  source  by  an  "equivalent"  source:  an  infinitely  thin  line  source,  covered 
by  the  foil,  and  having  a source  intensity  equal  to  the  average  effective  fraction 
of  the  gas  in  Fig.  3.5  times  the  actual  source  intensity.  This  fraction  we  find 

by  graphical  integration  to  be  f = 0.2116  at  R = 10  cm,  and  it  is  almost  inde- 

e 

pendent  of  R from  5 to  20  cm.  The  geometry  for  the  calculation  is  given  in 
Fig.  3.6. 

The  length  of  the  source  is  L(cm),  covered  by  a brass  foil  of  thickness 
t^(mil).  The  effective  source  intensity  is 


S = f S e/sec 
e e a 


emitted  isotropically  into  47T  steradians.  The  number  of  betas  emitted  per  sec- 
ond from  the  element  of  length  dX  is  the  energy  range  dE  (keV)  is 


j(E)dEdX  = (S  /L)|/(E)dEdX  (3.4) 

e 

where  j(E)e/(sec-keV-cm)  is  the  emitted  spectrum,  normalized  per  unit  length 
of  the  cylinder,  and  p(E),  which  contains  the  energy  dependence,  is  normalized 
to  unity.  From  Ref.  3.2  it  is  seen  that  the  spectrum  is  almost  flat  from  0 to  300  keV, 


"Cbscured  Point".  Located 
diagonally  opposite  a solid 
portion  of  cylinder. 


"Open  Point".  Located 
diagonally  opposite  portion 
of  cylinder  covered  only 
by  foil. 


Only  radiation  emitted  from  gas 
cylinder  in  shaded  regions  can  reach  "Point" 

in  question  without  scattering 


Fig.  3.  5 Geometry  for  Calculation  of  Effective  Gas  Fraction. 
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Fig.  3.6  Geometry  for  Analytical  Determination  of  Ionization 
Rate  from  Equivalent  Source. 


decreases  linearly  above  the  energy  E#  = 300  keV  and  is  zero  at  a maximum 
energy  E^  = 670  keV.  Thus  we  set 


|/(E)dE  = v dE 
o 


0<E<E 


= V dE(E  - E)/(E  -E  ) E $E<E 
o m m 8 8 m 


= 0 


E£E 


(3.5) 


m 


and  require  that 


rm 


*/(E)dE 
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which  yields 


completing  the  definition  of  V(E)  in  (3.5)  for  use  in  (3.4). 

The  betas  leave  the  element  dX  and  pass  through  the  foil  and  also  through 
a length  h(X)  of  air  of  density  P( g/cm3)  = 1.3xl0'3  g/cm3,  where 


Asa  result  the  energy  is  changed  from  E to  E\  and  the  flux  of  betas  emitted  by 
the  element  dX  in  the  energy  range  dE'  through  the  area  A'  at  the  point  of  inter- 
est is  F(E')e/(cm^-sec-keV-cm),  where 


47Th(X) 


The  energy  deposition  rate  E(keV-cm'  -sec'1)  at  the  point  of  interest,  due  to 
the  entire  beta  spectrum  and  source  length,  is 


S(E')t/(E)dE  (3.9) 


4rrL 


The  energy  integral  depends  upon  X and  t through  evaluation  of  the  stopping 
-12  1 

power  S(E')(keV-g  -cm  ) at  the  energy  E'(keV).  For  purposes  of  calculating 

z 

the  energy  loss,  the  path  length  traversed  is  expressed  in  g/cm  . Thus, 


R(g/cm  ) = PR(cm) 


Since  the  ratio  of  stopping  power  for  copper  (the  foil  is  actually  brass,  but  is 
essentially  copper  for  present  purposes)  to  that  of  air  is  almost  independent 
of  energy,  we  can  equate  the  foil  to  an  amount  of  air  that  would  produce  the 
same  energy  loss.  Use  of  S(E)  data  from  Ref.  3.3  for  Cu  and  air  yields 


, , c. 

tf(g/cm  air  equiv. ) = .01745  t (mils  Cu) 


(3.11) 


Hence,  the  path  length  T(g/cm  air)  for  a given  R is 


T = .01745  tf(mils)  + .0013R(cm) 


(3.12) 


and  the  path  length  equivalent  to  h(X)  in  g/cm  is  then  obtained  from  simple 
geometry  and  (3.12).  (See  (3.22)  below). 

The  stopping  power  of  air  (Ref.  3.3)  can  be  fitted  by  the  following  func- 
tion with  an  error  not  exceeding  about  +22%  from  10  to  900  keV, 


S(E')  = S (1  + E'/E  )' 
o o 


(3.13) 


where 


S = 5.778xl05  keV-g^-cm2 
o 

m = .5478 


E = .0135  keV 
o 


(3.14) 

(3.15) 

(3.16) 


This  stopping  power  expression  is  too  low  by  about  40%  at  1 keV,  correct  at 
0.1  keV  and  too  high  below  that.  However,  since  the  range  of  a 10  KeV  electron 
in  air  is  only  about  0.2  cm,  inaccuracies  in  S(E)  below  10  keV  are  unimportant 
for  the  present  work.  The  total  energy  of  each  particle  is  deposited  in  the  air 
surrounding  the  source  no  matter  what  S(E)  expression  is  used;  the  present  ex- 
pression could  lead  to  errors  of  0.2  cm  at  most  in  the  £(R,  v)  distribution. 

This  error  is  certainly  small  compared  to  that  due  to  neglect  of  scattering. 


18 


A significant  advantage  of  the  expression  (3.13)  is  that  it  can  be  inte- 

2 

grated  to  yield  the  distance  r(g/cm  air)  required  to  change  the  energy  from 
E to  E', 


r = f dE' 


/S(E') 


(3.17) 


[(l+E/Eo)m+1  -(l+E'/Eo)m+1]  (3.18) 

o 


The  range-energy  E^  is  that  value  of  E for  which  the  residual  energy  E'=0  at  r: 


Er(r)  = E0jf1  +So(m+l)r/Eo]1/(m+1)  - ij 


(3.19) 


Solving  (3.18)  for  E',  and  using  it  in  (3.13)  along  with  (3.19)  shows  that  S(E') 
can  be  written  in  terms  of  the  initial  energy  E and  the  range  energy  E 


S(E')  = S f(E) 
o 


(3.20) 


where 


f(E)  = [1  +(1  + E/E  )m+1  - (1  + E /E  )m+1]-m/(m+l)  (3  21 
o r o J 


E^  is  evaluated  for  that  path  length  r in  g/ciri  equivalent  to  the  distance  h(X)  in 
Fig.  3.6,  i.e. 


r (X,  v)  = T ■»  / 1 + (X-v)2/R2 


(3.22) 


where  T is  given  by  (3.12). 


- 3 - 1 

The  ion  -electron  pair  production  rate  Q(R,v)  ion  pairs-cm  -sec  is 


given  by 


Q(R,  v)  = £(R,  v)/I. 


(3.23) 


■where  T = .035  keV/ion  pair  is  the  energy  deposition  in  air  necessary  to  create 
one  ion-electron  pair.  By  use  of  (3.9)  and  (3.20)  this  result  can  be  written 
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1 


V 

where 


L 

Q(R,  v)/S  = / (Q  (R,v)/S  )dX 
a J x a 

0 


Q (R,  v)/S  = 
x a 


Pf  S 
e o 


4ttI  Lh(X) 


E 


m 

f(E)l/(E)dE 


(3.24) 


(3.25) 


The  lower  limit  on  the  energy  integral  has  been  replaced  by  E from  (3.19)  with 

r 

r given  by  (3.22),  since  no  particles  with  energy  less  than  this  can  reach  the 
point  of  interest.  The  result  (3.24)  is  called  the  specific  ionization  rate.  By 
converting  from  e/sec  to  curies  (Ci),  where 


1 Ci  = 3.7xl0^e/sec  (3.26) 

3 

the  units  of  Q/S  are  ion  pair s/ (cm  -sec-Ci),  and  that  of  Q /S  , the  specific 
a ^ x a 

ionization  rate  per  unit  length,  are  ion  pairs/(cm  -sec-Ci-cm). 

A FORTRAN  program  for  computing  the  result  (3.25)  by  double  numer- 
ical integration  is  given  in  Table  3.1.  The  program  input  and  typical  output  are 
given  in  Tables  3.2  and  3.3.  The  output  tabulates  (3.25)  and  (3.10)  vs  X and  pro- 
vides the  final  result  (3.24)  for  that  value  of  R following  the  table.  More  detailed 
results  are  given  in  Fig.  3.4,  labelled  "Analytical  Result  Neglecting  Scattering." 
In  general,  it  can  be  seen  that  the  calculations  agree  rather  well  with  the  ex- 
perimental data,  thus  reinforcing  the  general  accuracy  of  both  results. 

As  typical  values  for  use  in  the  charging  and  agglomeration  work  we 

choose  R = 1 0 cm  and  S = 500  Ci,  for  which 

a 

Q ~ 10^  ion  pairs/(cm^-sec)  (3.27) 


We  note  that  this  rate  drops  only  by  about  a factor  of  2 between  10  and  20  cm, 
so  that  a high  rate  is  maintained  throughout  the  volume  of  interest. 
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CYL ION . FT {FIND  ION  DEPOSITION  RATE  FROM  CYLINDRICAL 
BETA  RAY  SOURCE.  3/28/78. 

USE  AUG  OF  RATES  FROM  OPEN  AND  OBSCURED  SECTIONS. 

FRAC= .2116 
WRITE (1,108) 

READ  < 1 > 100 ) EL. » V , RMIN  , DR , RM AX 
R0=1 . 5095E- 08 
EMR=. 64608 
E0  = .0135 
READ ( 1 , 104)TF 
TF=1 . 745E-02*TF 

DUM=EO*<  ( 1 .+<TF+RMIN*1 . 3E-03 ) /RO ) **EMR~ 1 . ) 

WRITEC 1 , 105) DUM 
READ (1,101 )NX,NE,EMIN 
DX-EL/fL.OAT  ( NX ) 

EM1=1 .5478 
EMU=~. 35392 
ES=300. 

EM=A70 . 

EMS=370 . 

VL=1 . 303E+11/EL 
DE= ( EM- EM IN ) /FLOAT ( NE ) 

WRITE ( 1 » 107)DE 
R-RM IN-DR 
R=R+DR 

IF(R-RMAX)25,25» 10 
T-TF+R*1 .3E-03 
OEI-O 

WRITE(1» 106) 

DO  80  1 = 1 » NX 
X=DX*< FLOAT < I )- .5) 

SQ=SQRT(1 .+( <X- V)/R>**2> 

DUM=SG*T/RO 

ER=EO* ( ( 1 .+DUM)**EMR-1 . ) 

QE  J=0 

DO  60  J=1 > NE 

EJ=EMIN+DE*< FLOAT < J)~.5> 

IF ( EJ-ER >60,60, 30 

FJ=(  ( 1 .+EJ/EO)**EMl~LiUM)**EMU 

EN=1  . 

IF(EJ-ES)40»40,35 

EN=(EM-EJ)/EMS 

QEJ=QEJ+FJ*EN 

CONTINUE 

OEX=QE J*UL*DE*FRAC/ ( R*SQ ) **2 

WRITE  ( 1 ,102>X,GEX,ER 

QEI=QEI+QEX 

QEI=DX*GEI 

WRITE ( 1 » 103>R,GEI 

GO  TO  20 

FORMAT  < ' L= 'F5 . 1/ ' U= 7 F5. 1 / /RMIN= 'F5 . 1 / ' DELR= ' F5 . 1/ " RMAX= ' F5 . 1 ) 
FORMAT ( ,NX='I3//NE='I3/'EMIN=/F5.0) 

FORMAT (F5 . 1 » El  2 . 4 » F7 . 1 ) 

FORMAT ('  R='F5.1,'  CM  G/SA='E11.4'  I ONS/ ( CM**3-SEC-CURIE ) ' / ) 
FORMAT ( ' CU  FOIL  THICKNESS » MILS= ' F4 . 1 > 

FORMAT ( 'FOIL+MIN  AIR  PATH  RANGE  ENERGY , KEV= ' F6 . 1 ) 

FORMAT ( / ' X GX/SA  ER') 

FORMAT < ' DELTA  E='F4.1'  KEV'//'  X=CM, ' 

1'  OX/SA=IONS  / ( CM**3--SEC-CURIE-CM ) » ER-KEV'/) 

FORMAT <////' ALL  LENGTHS  ARE  IN  CM,  ENERGIES  IN  KFV ' / ) • 

END 


Table  3.  1 FORTRAN  Program  CYLION 
for  Calculation  of  Ionization  Rate  due  to 
Cylindrical  Beta  Source. 


21 


IRIS  PAGE  IS  BEST  QUALI'rv  PRiCTIfiAftbj 
FROM  COPY  FURNISHED  10  QDC  - 


R CYLION 


*> 


THIS  ?A9S  IS  BEST  QUJffilTI  TBX<SXW&* 
J8Q1I  oopx  C®SISHSD & W&  — ^ 


ALL  LENGTHS  ARE  IN  CM#  ENEPGIES  IN  KEV 

L = 30.  . 

V*1  5. 

PMIN*10. 

DELP-10. 

R MAX *20 . 

CU  FOIL  TH I CKNES  S.  M ILS* 2 . 0 

FOIL+MIN  AIR  PATH  RANGE  ENEFGY..KEV*  214.2 

NX*  30 

N E=2  30 

EMIN=2 1 0 • 

DELTA  E=  2.0  KEV 


X*CM,  QX/SA* I ONS/ < CM**  3- SEC-  CUFI E- CM ) t EP*KEV 


X 

QX/SA 

EP 

0.5 

0. 3301E+07 

308.8 

1 .5 

0. 3566E+07 

299.5 

2.5 

0.421 1E+07 

290.3 

3.5 

0.4680E+07 

281.2 

4.5 

0.5541E+07 

272.3 

5.5 

0 • 620  3E+07 

263.  7 

6.5 

0. 701 7E+07 

255.3 

7.5 

0. 7979E+07 

247.4 

8.5 

0.91 36E*07 

240.0 

9.5 

0. 1003E+08 

233.3 

10.5 

0.11 1 0E+08 

227.3 

11.5 

0. 1 242E+08 

222.  3 

12.5 

0.1 333E+08 

218.4 

13.5 

0. 1 368E+08 

215.7 

14.5 

0. 1 428E+08 

214.4 

15.5 

0. 1428E+08 

214.4 

16.5 

0.1 368E+08 

215.7 

17.5 

0.1 333E+08 

218.4 

18.5 

0.1 242E+  08 

222.3 

19.5 

0.1 1 10E+08 

227.  3 

20.5 

0. 1003E+08 

233.  3 

21 .5 

0.91 36E+07 

240.0 

22.5 

0. 7979E+07 

247.4 

23.5 

0.701 7E+ 0 7 

255.  3 

24.5 

0.6203E+07 

263.  7 

25.5 

0 • 554 1 E+  07 

272.  3 

26.5 

0 . 4680E*  07 

281.2 

27.5 

0.421 1 E+07 

290.3 

28.5 

0 . 3566E+07 

299.5 

29.5 

0. 3301E+07 

308.8 

Table  3.  2 Input  to  Program  CYLION  and 
Output  for  R = 10  cm. 


R*  10.0  CM  Q/SA*  0 . 2530E+09  IONS/ C CM** 3- SEC- CUP IE) 
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1 


i 4# 


♦J* 

* # 


«SSSSE*»,,# 


X 

QX/SA 

ER 

0.5 

0.1 844E+  0 7 

286.7 

1.5- 

0.  1 927E+07 

262.4 

2.5 

0. 2040E+07 

278.3 

3.5 

0.21 72E+07 

274.3 

4.5 

0»  2 334E+07 

270.6 

5.5 

0.2346E+07 

267. 1 

6.5 

0.2515E+07 

263.9 

7.5 

0.3147E+07 

261  .0 

8.5 

0.2788E+07 

258.4 

9.5 

0.2859E+07 

256. 1 

10.5 

0.29S4E+07 

254.2 

11.5 

0 • 3087E+0  7 

252.6 

12.5 

0. 3024E+07 

251.4 

1 3.5 

0. 3187E+07 

250.6 

14.5 

0.31 49E+07 

250.2 

15.5 

0.3149E+07 

250.2 

16.5 

0*316  7E+  0 7 

250.6 

17.5 

0. 3024E+07 

251.4 

18.5 

0. 308  7E+07 

252.6 

19.5 

0 • 2954E+07 

254. 2 

20.5 

0 • 2859E+  0 7 

256.  1 

21.5 

0 • 2788E+07 

258.4 

22.5 

0.3147E+07 

261.0 

23.5 

0.251 5E+07 

263.9 

24.5 

0.2346E+07 

267.  1 

25.5 

0 • 2 334E+  0 7 

270.6 

26.5 

0.21 72E+07 

274.  3 

27.5 

0 . 2040E+  0 7 

278.3 

28.5 

0. 1927E+07 

282.4 

29.5 

0.1 844E+07 

286.7 

Table  3.  3 Output  of  Program  CYLION 
for  F = 20  cm. 
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4.  CHARGE  GENERATION  AND  SEPARATION 


The  generation  and  separation  of  charge  takes  place  in  one  cylindrical 

section,  as  shown  in  Figs.  2.1  and  2.2.  The  idealized  geometry  of  one  segment 

85 

of  this  cylinder  is  shown  in  Fig.  4.1.  The  Kr  produces  ionization  in  the  cen- 
tral portion  of  the  segment,  with  a potential  (volts)  across  the  section  edge 
plates.  For  simplification  in  the  theoretical  analysis  of  the  charge  separation, 
the  reduced  variable  y (range  0 to  1 ) is  introduced,  and  plane  parallel  geometry 
is  assumed  with  the  separation  d - R0  being  a (variable)  parameter. 


Fig.  4.1.  Geometry  of  Charge  Generation  and  Separation 
Section. 


The  charge  density  and  electric  field  equations  can  be  written  as 


\aod2/  dy2  \Qod/  dy  Q0  \ Qq  / + ' 

(4nepn  2\ 

-Q^-)  N+(N+-  NJ 


Ci  \ d^N. 


^EnQ\  dN_ 


Qod2/  dy2  \Qod  / dy 


Q fanQZ\ 

~ - — ^ N N 

Ri  \Qo  / + ' 


4u  ep  n 


N_(N+  - NJ 


= (4n  en0)(N+  - N J 


where  D is  the  diffusion  coefficient,  — 1 0*V(cm3- sec)  is  the  estimated  ion 

— — ft  1 

pair  production  rate  at  (typically)  R = 10  cm  (Sec.  3),  n = \/Q/oc  w 2. 4x10  /cm3  is 

o * o 

the  zero  field  equilibrium  ion  density  for  a recombination  coefficient 

•63  2. 

cr=*1.7xl0  cm  /sec  for  air,  p » l cmV(V-sec)  is  the  ion  mobility  in  air,  and 
N+  =n+/nQ  are  the  normalized  positiveAiegative  ion  densities,  E is  the  elec- 
tric field,  and  d and  y are  as  defined  in  Fig.  4.1.  We  define 

pn  V ftV  1 

, o o _ 'o 


° d*/o~d 


=*  V^(voltsV41 231 , for  d = 10  cm,, 


neglect  the  diffusion  term  (D/d«pE  for  reasonable  E),  and  obtain 
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n 


»* 


- » 


dN 

E + Q 

0 — ~ - N N 

o E dy  Q + - 
o 7 o 


E_  dN-  Q 

E dy  * Q N+N- 
o o 


PN+  + PN+N_ 


- PN 


+ PN+N_ 


(4.5) 


(4.6) 


where  P = 4«pen  /Q  = 4 it  p e/a  =»  1.06  and  E = V /d. 

o o o o 

The  N+N_  terms  in  (4.5)  and  (4.6)  drop  out  if  P = 1 , and  since  P is  very 
nearly  1.0,  we  set  P - 1.0  and  obtain  a simplified  set  of  equations  (an  approach 
consistent  with  the  present  feasibility  determination): 


o E dy 
o ' 


dN 


Q 

Q 


o E dy  Q 
o 7 o 


N 


N 


and 


d2^ 

72  = N-  - N+ 

dy 


(4.7) 


(4.8) 


(4.9) 


where  pn^V/fQ^d  )>  with  V the  potential  at  position  y. 

Consider  the  situation  of  Fig.  4.1,  where  0=0  for  Osy<y^r  and  y^ysl, 
and  Q=Qq  for  yjsysy^  For  this  case  the  solution  to  (4.7),  (4.8)  and  (4.9)  can 
be  written  as: 


0- y-Yj : 


N+  = 0 


N-  = w *f  z 


d^_ 


O 0 


d^/dy 


(4.10) 


ylsysy2: 


N+  = tanhf  f — 
^1 


Jx. 


d^/dy 


(4-11) 


N-=tanh[/ 
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y2-y-1: 


1 = [“  + fY 

+ N J d^/dy 

*■  ° V-J 


N = 0 


(4.12) 


where 


N = tanh 
o 


[A 


- fo_  . 

d^/dy 


(4.13) 


The  set  of  equations  (4.10)  to  (4.13)  can  be  used  to  solve  for  N+I  N_,  and 
all  as  functions  of  y,  by  splitting  O^ysl  into  several  segments  and  numerically 
solving  for  the  values  at  each  breakpoint.  This  has  been  done  for  16  break- 
points (including  y=0  and  1 ),  and  the  results  for  N+  and  N_  are  plotted  in  Fig. 
4.2  for  three  values  of  <t>Q.  Only  one  N_  profile  has  been  plotted,  since  all 
N_  profiles  are  the  N+  profiles  reflected  about  y = 0. 5.  Note  that  the  results 
imply  that  some  optimum  value  of  exists  which  maximizes  N+(y=l). 

The  net  charge  separation  (N+  - N ) is  plotted  in  Fig.  4.3.  As  a result 
of  the  method  of  normalizing  to  nQ,  the  results  do  not  depend  explicitly  onQ, 
the  ion- electron  production  rate.  Q must  be  determined  analytically  and  ex- 
perimentally to  allow  more  accurate  calculations  for  N+,  N . Figure  4.3 
shows  that  while  small  iQ  maximizes  the  charge  separation  at  y = y ,y  , the 

X u 

maximum  integrated  charge  separation  is  for  some  intermediate  value  of  6 . 

o 

The  best  value  of  4 is  that  which  maximizes 
o 


N = f (N  - N )dy 
sep  J + - 7 


(4.14) 


and  this  can  be  readily  done  by  numerical  calculations.  Because  of  the  nor- 
malization of  N and  N to  n , ^here  is  a single  optimum  for  ^ , and  the  cor- 

T — O O 

responding  can  be  readily  determined  from  Eq.  (4.4)  for  any  set  of  para- 
meters d,  etc.  Since  d and  n^  = ^Q^/a  vary  somewhat  with  R (see  Fig.  4.1), 
the  optimum  varies  slightly  with  R,  and  a suitable  average  value  will  have 
to  be  selected.  From  Fig.  4.3  it  can  be  seen  that  net  charge  concentrations 


v 


of  >0.5  nQ  can  be  achieved  in  the  separation  regions,  so  that  separated  charge 
densities  of  more  than  half  the  equilibrium  electron-ion  density  n^  can  be 
achieved.  That  these  densities  will  be  sufficient  to  induce  an  amount  of  charge 
on  particulates  that  will  cause  significant  agglomeration  in  that  section  of  the 
device  is  investigated  in  Section  5. 

The  geometry  of  Fig.  4.1  is  idealized,  since  in  practice  the  ionization 
region  will  not  have  sharp  boundaries.  Indeed,  the  results  of  Section  3 show 
that  beta  particle  scattering  produces  very  significant  effects.  However,  the 
theoretical  analysis  shows  that  significant  charge  separation  can  be  achieved 
and  gives  numerical  values  for  the  approximate  optimum  potential  Vq.  This 
helps  significantly  in  guiding  any  experimental  work  done  in  the  next  stage  of 
work. 


v 
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5.  PARTICLE  CHARGING  AND  AGGLOMERATION 

In  this  section  we  consider  how  the  particles  are  charged  in  the  charger 
(Sec.  5.1)  and  how  they  are  then  agglomerated  in  the  agglomerator  (Sec.  5.2). 

We  assume  (a)  that  there  are  essentially  two  populations:  "small1'  particles 
below  about  one  micrometer  in  radius  (subscripts  i ),  and  "large"  particles 
above  about  one  micrometer  in  radius  (subscripts  2);  and  (b)  that  the  small 
particles  agglomerate  onto  the  large  particles. 

5.1  Charging  Section 

In  the  charging  section  we  assume  that  the  particles  find  themselves  in  a 
charging  environment  where  there  are  unipolar  ions  of  density  n+,  and  an  elec- 
tric field  of  intensity  E^.  We  assume  the  particles  are  exposed  to  this  envir- 
onment for  an  interval  of  one  second.  (The  results  are  insensitive  to  the  exact 
value  of  this  time  interval.)  For  small  particles,  the  charging  is  generally 
dominated  by  the  "diffusion"  mechanism,  while  for  large  particles  the  charging 
is  generally  dominated  by  the  "field"  (conduction)  charging  mechanism  (Ref.  5. 1,5. 2). 
For  the  present  purposes,  we  assume  here  that  the  small  and  large  particles 
are  charged  only  by  diffusion  and  conduction,  respectively. 

5.1.1  Diffusion  Charging 

For  diffusion  charging,  which  is  appropriate  for  small  particles,  the  num- 
ber of  elementary  charges  acquired  in  time  t is  given  [Ref.  5.1,  Eq.  (3.61); 

Ref.  5.2,  Eq.  (5.27)]  by: 


RkT 

2 


1"(‘ + k) 


wher  e 


kT 


7Te*U.n  R 
1 + 


(5.1) 


(5.2) 


where  R is  the  radius  of  the  particle,  n+  is  the  undisturbed  ion  density,  e is  the 
magnitude  of  the  electron  charge,  U.  is  the  rms  thermal  velocity  of  the  ions, 
k is  Boltzmann's  constant,  and  T is  the  absolute  temperature.  At  room  tem- 
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perature,  we  take  kT  = 4x10  erg  and  IR  = 5x10  cm/sec  for  nitrogen  mol- 
ecules, and  we  may  write: 

_ef_  _ 5.76xl0~6 

RkT  R '5 


TD  0.905n+R  8eC 


n,  = 1.736x10  Rln(l  + 0.905n+Rt) 


(5.5) 


Assuming  (a)  the  optimum  value  of  for  separating  the  positive  and  negative 
ions  (given  by  <f>^  - 0.378  in  Section  4),  (b)  a distance  of  d = 10  cm  = arc  length 

between  the  plates,  and  (c)  a local  source  ionization  rate  of  Q = lO^ion  pairs 

, * 8 - 3 ° 

per  cmJ  per  sec,  we  obtain  a value  n sr  1.5x10  cm  for  ion  density,  at  a volt- 
age of  about  15600  volts  across  the  plates  and  a field  intensity  of  E = 5.2  stat- 
volts/cm.  (The  value  of  n+  is  the  average  value  of  the  ion  density  profile,  that 
is,  about  60%  of  the  nominal  ion  density  n^  = y/Q^/Of,  where  Qq  is  the  ionization 

rate  = lO^cm  ^sec  *,  and  a is  the  recombination  coefficient  = 1 . 7x1  0 ’^cm^iec  ^.) 

8 - 3 

For  particles  of  radius  R = 0.1  micrometer,  and  n+  = 1.5x10  cm  , Eq.  (5.1) 
becomes 

n,  = 1.736  ln(l  + 1 358t)  (5.6) 

= 12.52  charges  in  one  second. 

Note  that  this  charging  is  insensitive  to  the  time  t and  to  the  ion  density  n^  be- 
cause of  the  logarithm.  It  is  explicitly  independent  of  the  charging  field  E 
but  n+  depends  on  this  as  seen  in  the  previous  section. 


5.1.2  Conduction  (or  Field)  Charginc 


For  conduction  (or  "field")  charging,  which  is  appropriate  for  large  par- 
ticles, the  number  of  elementary  charges  of  magnitude  e acquired  in  time  t is 
given  [Ref.  5.1 , Eq.  (3.  58);  R ef.  5.2.  Eq.  (5.15)]  by: 


' nsat(t  + T ) 
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«> 


wher  e 


and 


r = — - — = 2.210xl06/n  sec 
c TTe^n  + 


n = 3R2  E /e  = 6.25x10  E , R2 
sat  ch  ch 


(5.8) 


(5.9) 


where  n denotes  the  "saturation"  charge,  a finite  value  attained  after  an  in- 
finite time,  with  the  time-scale  given  by  r . Equation  (5.8)  for  T is  obtained 

c c 

assuming  = 1 cm2/volt-sec  = 300  cm2/sv-sec.  For  practical  purposes,  n = nBaj 

may  be  considered  the  value  of  n2  after  a few  time  periods  t . Note  that  n does 

c sat 

not  depend  on  n+,  while  T essentially  depends  only  on  n and  is  independent  of  Rand 

8 3 

Ech*  For  ion  density  n+  = 1.5x10  /cm  , Tc  is  only  0.015  sec  so  that  on  a time 

scale  of  the  order  of  seconds  all  large  particles  have  n2  s n . (The  factor  3 

sat 

is  Eq.  (5.9)  is  for  conducting  particles;  for  nonconducting  particles  this  factor 
is  between  unity  and  3,  depending  on  the  dielectric  constant.)  For  R = 0.1 //m 
and  Ech  = 5*2  sv/cm,  n2  = 3.25  by  conduction  compared  with  n,  -12.52  obtained 


with  one  second  of  diffusion  charging.  For  R = 1.0/;m  and  the  same  E , n2  =325 

ch 

by  conduction  char ging  compared  with  n,  =165  by  diffusion  char  ging.  At  R=0.468 f^xn 


we  have  n,  = n2  = 71.1.  Hence,  the  transition  between  diffusion  and  conduction 
as  the  dominant  charging  mechanism  occurs  in  the  range  R = 0.  l//m  to  R = l.Ofxm. 
Thus,  both  mechanisms  are  important  in  contributing  to  the  charge  in  this  crit- 
ical range.  However,  since  we  are  interested  in  all  values  of  R below  l.O^m 
(all  submicron  particles),  we  will  be  on  conservative  ground  if  we  assume  that 
diffusion  alone  is  responsible  for  the  charging  of  all  submicron  particles. 

The  charging  theories  used  for  the  estimates  of  this  report  are  the  standard 
simplified  theories.  There  exist  more  complicated  theories  which  combine  the 
diffusion  and  conduction  mechanisms  (Ref.  5.3),  but  *heir  use  would  not  change 
the  numbers  significantly  and  seems  unjustified  for  our  present  purposes. 

The  large  particles  are  assumed  to  act  as  collectors  and  to  be  described 
by  a log-normal  distribution.  Among  the  quantities  needed  in  the  theory  of  the 
*8f?lorneration  process  to  be  considered  next  are  the  average  values  of  n and  of 
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r 


«n» 

■*  / 


r. 


Bn,  where  B is  the  particle  mobility  assumed  to  be  given  by  the  Stokes  viscous 
drag  law: 


B = C[R)  (*££.) 

6nr)R  \ g / 


(5.10) 


»-4 


where  r)  is  the  dynamic  viscosity  of  air  (which  we  take  to  be  1.8x10"  poise  or 
g/cm/s ec  at  room  temperature),  and  C(R)  is  the  Cunningham  "slip"  factor  which 
accounts  for  the  fact  that  the  particle  is  not  large  compared  with  the  collisional 
mean  free  paths  of  the  air  molecules.  The  latter  factor  has  been  given  as 
(Refs.  5.4,  5.5): 


C(R ) = 1 + [1.257  + 0.4  exp  (-  1.1R* )]  /R ' 
with  R'  = R/6.62xl  0 ^ = R(//m)/  0. 0662 


(5.11) 


Table  5.1  shows  how  C(R)  varies  with  R.  Above  R = 1.0/im,  C(R)  is  essentially 
unity.  Below  R = 0.1/jm,  C(R)  rises  approximately  as  R_*.  Table  5.1  also 
shows  the  product  of  B and  n (Bjn,  for  "small"  particles),  where  n is  obtained 

Q 2 

from  the  diffusion-charging  formula  [Eq.  (5.5)]  with  n+  = 1.5x10  /cm  and 
t = 1 sec.  It  is  important  to  note  that  the  product  Bn  is  almost  constant  in  the 
range  R = O.l^m  to  R = 1.0//m.  Thus,  since  the  agglomeration  theory  to  be 
described  next  involves  n,  through  the  product  B,  n,  for  small  particles,  the 
results  will  be  insensitive  to  R,  in  the  range  0.1^/m  to  1.0|/m. 

5.2  Agglomerator  Section 

The  air  flowing  into  the  agglomerator  section  is  loaded  inhomogeneously 
with  charged  particles,  such  that  there  are  essentially  only  positive  charges 
on  one  side,  and  negative  charges  on  the  other  side.  There  is  a superimposed 
electric  field  of  intensity  E in  this  section  whose  purpose  is  to  pull  the  oppos- 
itely-charged clouds  of  particles  toward  each  other  and  to  mix  them,  so  that 
oppositely-charged  partners  can  come  together  and  agglomerate.  The  task  of 
this  section  is  to  determine  the  time-scales  under  which  the  two  processes  occur: 
transit  across  the  gap,  and  agglomeration. 
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TABLE  5.  1 


Variations  of  C(R)  and  B,nj  with  R* 

8 ^ 

(Assumes  n+  = 1.5x10  /cm  , t = 1 sec  diffusion  charging) 


R (pm) 

C(R ) 

B n 
l l 

B 

l 

n 

l 

n 

sat 

0.0 

• 

00 

7.620x10  10 

00 

0 

0 , 

0.001 

110. 26 

1.5 12x10 10 

3.250xl0U 

0.0465 

0.000325 

0.01 

14.93 

3. 757xl09 

4.400xl09 

0.854 

0.0325 

0. 1 

1.882 

6.949x10^ 

5.547xl0? 

12.52 

3.25 

0.2 

1.421 

5.749xl08 

2. 094x10  7 

27.45 

13.0 

0.3 

1.278 

5.436xl08 

1.256xl0? 

43.29 

29.25 

0.4 

1.208 

5.3l6xl08 

8. 90  lxlO8 

59.72 

52.0 

0.5 

1. 166 

5.266xl08 

6.873xl06 

76.58 

81.25 

0.6 

1.139 

5.247xl08 

5.595xl06 

93.80 

117.0 

0.65 

1.128 

5.244xl08 

5. 115xl06 

102.5 

137.3 

0.7 

1.119 

5.244xl08 

4.7  V lxlO8 

111.3 

159.3 

0.8 

1.104 

5.250xl08 

4.067x10^ 

129.1 

208.0 

0.9 

1.093 

5.26  lxlO8 

3.579xl06 

147.0 

263.3 

1.0 

1.083 

5.275xl08 

3.192xl06 

165.2 

325.0 

2.  0 

1.042 

5.44 lxlO8 

1.536xl06 

354.5 

1300. 

*R  = particle  radius 
C(R)  = Cunningham  slip  factor 
B j = particle  mobility 

n,  = particle  charge  (number  of  elementary  charges)  by  diffusion 

n = saturation  charge  by  conduction  at  E = 5.2  sv/cm 
sat 


v 

J 
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5.2.1  Agglomeration 

We  wish  to  consider  two  clouds  of  oppositely-charged  particles,  moving 
past  each  other  with  constant  velocity,  andto  determine  the  probability  of  a con- 
tact occurring.  In  particular,  we  consider  a large  particle  acting  as  a collec- 
tor, with  a reference  frame  fixed  on  this  collector.  The  other  particles  stream 
past  it  at  velocity  v^,  so  that  we  can  determine  the  collection  efficiency  by  analyz- 
ing the  orbits  o‘f  the  streaming  particles  and  counting  those  which  intersect  the 
collector.  To  do  this  analytically,  we  assume  that  the  collected  particles  are 
inertia-less,  that  is,  their  inertial  acceleration  is  negligible  compared  with, 
say,  that  due  to  viscous  drag  (demonstrated  by  calculations  in  Ref.  5.5;  see 
also  Ref.  5.2,  p.  36).  An  appropriate  analysis  for  this  problem  approximates 
the  orbit  by  that  resulting  from  a constant  velocity  v^  superimposed  onthe  velocity 
due  to  a central  force.  If  these  velocities  are  sufficiently  large,  it  can  be 
shown  that  Brownian  diffusion  can  be  neglected  (Ref.  5.5).  (That  is,  including 
Brownian  diffusion  will  lead  to  the  same  result  in  the  limit  of  large  charges  - 
see  Sec.  5. 1 . 2. ) 

In  Fig.  5.1,  we  depict  the  collection  of  small  particles  of  radius  R,  by  a 
large  particle  of  radius  R2  , under  an  attractive  central  force  of  magnitude 
Kp/r*3  between  the  particles,  with  p denoting  the  exponent  for  a power-law 
force,  and  Kp  the  associated  constant;  v^  denotes  the  relative  velocity  = 

(B^j  + B2  n2  )eE.  There  are  3 types  of  orbits  portrayed,  Type  A represent- 
ing uncollected  particles,  Type  B representing  collected  particles,  and  Type 
C representing  the  critical  orbit  dividing  the  two  classes  of  orbits.  If  we  can 

determine  x , the  initial  offset  distance  of  the  critical  orbit,  then  the  collec- 
o 

tion  efficiency  is  given  by  E-  x 2/(R.  + R2  )2.  To  compute  x we  write  the 

o 1 o 

equations  for  the  orbits  in  the  form 


v = v-  - (K  /rP)  r 
o p 


(5.12) 


~ = v = -v  cobO  - K /rP 
dt  r o p 


(5.13) 


'-tv  , A 

rTt  = vo  =vo8in* 


(5.14) 


From  these  equations  the  shape  of  the  orbit  is  given  by 


ldr  . . „ *'p 

- — - - cote  

r d0  p . . 

v r sin© 
o 


(5.15) 


from  which  we  obtain  the  solution 


pK 

P p 
x = x - 

O V 


/0 

5inp-S 


(5.16) 


where  x = rsine 

The  critical  orbit  is  that  for  which  x vanishes  when  0 = 7T , that  is,  the  right- 

hand  side  vanishes.  Thus,  we  have  the  critical  value  of  x : 

o 


T r” 

, . „p->. 

o v f 

° J 

0 


(5.17) 


and  the  collection  efficiency  may  be  written: 


x 2 pK 

c 2 \ _ J 

Cp  (R,+R,)2  ~ (R,  +R2)2  vq 


K f*  2 

HWsinP-W 

0 


(5.18) 


For  a Coulomb  force  of  attraction  between  the  two  particles  (the  case  of  inter- 
est here),  p and  are  given  by  p = 2 and  K^=B,  n,n2  e2  , where  n,  and  n2  are 
the  numbers  of  elementary  charges  on  the  particles.  The  collection  efficiency 


4B,  n,  n2  e2 
? = vo(R1+R,)2 


(5.19) 


If  and  7} i denote  the  density  of  small  particles  and  large  particles,  res- 
pectively, the  rate  of  disappearance  of  small  particles  per  unit  volume, dT^j/dt, 
is  given  by  multiplying  by  and  by  the  current  of  small  particles 
*(R,+R2)2  v 


v 


^ = - 4 it%  <7}l  Bj  n,n2  e4  = - 


(5.20) 


Thus,  rA  represents  the  agglomeration  time 


TA  = (4^2B,n,n2.2)_1 


(5.21) 


which  is  seen  to  be  independent  of  E (through  vq).  Assuming  7]  z constant, 
this  implies  an  exponential  decay  of  the  small -particle  population,  with  time- 
constant  T , due  to  agglomeration  onto  the  large  particles. 

If  the  small  particles  have  no  charge,  they  are  still  attracted  to  the  charged 
large  particles  by  the  polarization  force  which  is  inversely  proportional  to  the 
fifth  power  of  the  inter-particle  distance.  (R  ef.  5.5).  In  this  case  we  set  p=5  in 
Eq.  (5.18),  and  K^  = 2Rj  Bjn2  ^ e4.  Evaluating  the  integral  and  deriving  an  equa- 
tion of  form  similar  to  Eq.  (5.20),  we  find  that  the  agglomeration  time  for  the 
neutral-small-charged-large-particle  case  is  given  by 


■ft  - 1 1 5 “v 

= 7*  0 


7/2  3/2 


B,R,  n2V 


(5.22) 


which  is  (of  course)  independent  of  n,,  and  which  depends  on  E through 
V (B»ni  + Bjn*)eE.  This  2/5  -power-law  result  corresponds  to  the  asymptotic 
limit  of  results  obtained  previously  by  Parker  (Ref.  5.5)  using  boundary-layer 
analysis  to  take  into  account  fluid-flow  effects  in  the  collection  of  uncharged 
aerosols  by  charged  droplets. 

5.2.2  Thermal  Agglomeration 

The  formula  for  agglomeration  time  T^,  Eq.  (5.21),  was  derived  by  assum- 
ing inertia-less  particles  and  neglecting  diffusion.  Zebel  (Ref.  5.6)  treats  the  ag- 
glomeration of  charged  aerosol  particles  undergoing  Brownian  diffusion.  For 
two  types  of  particles,  of  radius,  mobility,  and  charge  R,,  B,,  and  n,e.  and 
Rit  B2,  and  n2  e,  the  ratio  of  thermal  agglomeration  rates  (or  "coagulation 
constants")  for  charged  and  uncharged  particles  is  a factor 
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(5.23) 

(5.24) 

and  Go  is  a 8as;- kinetic  correction  factor  which  vanishes  for  large  particles  and 
which  introduces  finite  mean-fr ee-path  effects  for  small  particles.  The  variable 
y describes  the  electrical  force  between  the  charged  particles;  it  vanishes  for 
neutral  particles  and  is  large  for  large  forces,  positive  for  repulsion  and  neg- 
ative for  attraction.  For  strongly-charged  aerosols  (|y|  large),  the  increase 
due  to  attraction  exceeds  the  decrease  due  to  repulsion  and  results  in  a net  sig- 
nificant increase  in  agglomeration.  For  large  negative  y,  f becomes  equal  to  |yj. 
Thus,  since  the  neutral  agglomeration  rate  per  unit  volume  is 

Ko  = 47T(B,  + B2  )(R,  + R2  )kT?7,  7]l  (5.25) 

the  net  agglomeration  rate  is 

Kof  = 47t(D,  + B2  )n,n2e^  ?/,  7jt  (5.26) 

which  corresponds  to  the  right-hand  side  of  Eq.  (5.20)  if  we  neglect  B2  compared 
with  Bj  (which  should  be  justified  if  B2  and  B]  are  associated  with  large  and 
small  particles,  respectively). 

5.2.3  Electrostatic  Dispersion  or  Mutual  Repulsion 

Another  effect  which  needs  to  be  considered  is  that  due  to  the  charges  within 
the  two  oppositely-charged  clouds  of  unipolar  particles  before  they  are  brought 
together  in  the  agglomerator . Each  cloud  will  tend  to  expand  and  disperse  to  the 
walls  as  the  particles  recede  from  each  other  under  their  mutual  repulsions.  In 
particular,  we  desire  to  know  the  time-scale  under  which  the  large  particles, 


which  are  to  act  as  collectors,  will  disperse.  The  decrease  in  concentration  of 


4 4 

o 


where 


f = 


(1  + Goy)exp(y)  -1 


n,  n2  e 


(R  i +R2)kT 


•*** 


a spatially  - homo  g en  eou  s cloud  of  large  particles  (but  with  distributed  sizes) 
is  given  by 

~ - div(^/2  v)  = - 71 2B2  n2  e divE 

= -‘V2B2r\2  e-47TeZ^?n 
j j J 

= -?/2B2n2  ■ 4ne2  7]  n 
o 

= -^/tr  (5.27) 

where  7]  is  the  total  number  of  particles  per  unit  volume  and  n is  the  average 
number  of  charges  on  a particle.  Thus,  the  average  repulsion  time  of  a cloud 
of  distributed  sizes  and  charges  may  be  taken  to  be: 

Tr  = (47T??oB2n2 -He2)'1  (5.28) 

where  B2n2  is  the  average  value  of  B2n2  in  the  cloud. 

Thus,  the  particle  clouds  must  be  brought  together  within  a time  less  than 

TR  *n  ort*er  to  8et  them  mixed  so  they  can  agglomerate.  Once  mixed,  they  will 

require  a time  of  order  T to  agglomerate  significantly.  Formulas  for  T and 

A A 

Tr  useful  for  system  design  purposes  will  be  presented  in  the  next  section, 

where  particle  size  distributions  are  taken  into  account. 

5.3  Particle  Size  Distributions  and  Design  Criteria 

Although  theoretical  analyses  are  simpler  when  the  particle  size  distribu- 
tions are  monodisper se,  in  actual  cases  it  is  necessary  to  deal  with  polydis- 
perse  aerosol  systems.  Such  systems  are  usually  described  by  log-normal 
distributions.  The  properties  of  log-normal  distributions  are  well  known  (e.g. , 
Ref.  5.1,  Chap.  5).  In  this  section  we  will  replace  functions  of  particle  diameter 
D by  mean  values  of  these  functions  based  on  the  log-normal  distribution. 
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V 


Typically,  particle  distributions  are  defined  by  the  mass  median  diameter 
(here  denoted  by  D^^)  and  the  geometric  standard  deviation  (here  denoted  by 
Og).  All  other  special  quantities  of  interest  may  be  derived  from  these.  Among 
the  special  quantities  are  the  following: 

D = count  median  or  geometric  mean  diameter  = D exp  (-3  In2  a ) 
g m /2  g 

D = count  mean  diameter  (or  simply  "mean  diameter")=  D exp(0.51n2a  ) 

g g 

— = count  mean  square  diameter  (average  area)  = D 2 exp  (2  In2  a ) 

g g 

— 5 

D = count  mean  cube  diameter  (average  volume)=  D ^exp  (4.5  In2#  ) 

g g 

F , = relative  number  (or  fraction)  of  particles  with  diameter  D less  than  D, 

(fraction  of  "small"  particles) 

= 0.5  erfc  [(In  D - In  D,)/y21n  a 1 

g r 

Fz  = relative  number  (or  fraction  of  particles  with  diameter  D greater  than  D, 
(fraction  of  "large"  particles) 

= 0.5  erfc  [(In  D,  - In  D )//7ln  a ] 

g g 

= 1 - F, 


5.3.1  Agglomeration  Time 

The  agglomeration  time  is  given  by  Eqs.  (5.20)  and  (5.21).  We  will  re- 
place 7}  2 by  F where  7/  is  the  total  number  of  particles  per  cm^,  and  7] 


n - W lM- 

° f pP  * 


(5.29) 


where  77]  denotes  the  mass  loading  in  grams  of  aerosol  per  cm^  of  air,  p is  the 

solid  density  of  the  aerosol  (here  taken  to  be  2 gm/cm^),  and  is  the  mean 

cube  of  the  particle  diameter.  We  also  replace  n,  by  the  mean  value  of  n , 

sat 

which  is  (Eq.  (5.9)): 


n = 3R*E  . /e  = L/4e 

sat  ch  ch 


(5.30) 


where  D2  is  the  mean  square  of  the  particle  diameter.  Then 


may  be  written: 


(5.31) 


where  the  second  parentheses  depends  only  on  the  distribution.  Note  that  B,n, 
depends  only  on  D,  (=  2Rj),  and  that  T is  inversely  proportional  to  the  mass- 

A 

loading  parameter  We  take  P-  2 gm/cm3  and  E , = 5.2  sv/cm,  and  assume 

8 Ch 

R - 0. l/zm  so  that  B , n,  = 6.949x10  (Table  5.1).  Then  Eq.  (5. 31 ) may  be  written: 

rA^.6.406xl0-6(?^sr)  (5.31) 

(R  = 0.ty/m)  in  fim 

It  is  important  to  note  that,  while  T is  evaluated  here  only  for  R = 0. 1/lm,  the 

A 

dependence  on  R is  through  Bjnt  , which,  as  Table  5.1  shows,  is  nearly  constant 
in  the  R -range  O.lfi m to  l.Ofpn. 


5.3.2  Repulsion  Time 

The  mutual  repulsion  time  T is  given  by  Eqs.  (5.27)  and  (5.28).  We  re- 

K. 

place  ?7o  by  Eq.  (5.29),  and  n by  n as  in  Eq.  (5.30).  The  average  value  of 
B2n2  is  replaced  by 


B 


- - (-L—  3EchD^ 
\3nr]D  ' 4e  / 


ch 

4neri 


Then  may  be  written: 


= 9.295x10 


-6 


(5.33) 


(5.34) 


where  we  have  assumed  P=2,  ^=1.8x10“^,  and  E , =5.2,  and  where  the  final 

ch 

parentheses  depends  only  on  the  distribution.  As  in  the  case  of  T , T is 

A R 
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inversely  proportional  to  % The  dependence  on  inverse  E * results  from 

ch 

the  assumption  that  for  the  purposes  of  estimating  repulsion  times  all  particles 
charge  up  as  large  particles,  i.e.,  by  conduction  alone.  This  should  yield  a 
good  estimate  of  since  most  of  the  charge  resides  on  the  large  particles, 
and  the  error  in  assigning  the  wrong  value  of  charge  to  the  small  particles 
should  be  small. 

Note  that,  after  the  mixing  has  been  achieved,  oppositely-charged  large 
particles  will  tend  to  agglomerate  with  one  another.  This  large-particle  ag- 
glomeration time  is  roughly  given  by  T . 

R 

5.3.3  Voltage  to  Bring  Oppositely  Charged  Clouds  Together 

The  voltage  required  to  force  the  oppositely-charged  clouds  together  a- 
crossa  gap  of  length  d may  be  written  as: 


V(volts)  = 


300  d2 

(B  i n j + B2  n2  )et 


300  d* 

[B,n,  + DE  ,/(4neT])]ct 


(5.35) 


where  we  have  used  Eq.  (5.3  3)  for  B2n2,  and  the  factor  300  is  required  to  con- 
vert statvolts  to  volts.  The  time  t should  be  T orT  , whichever  is  shorter. 

A R 

We  assume  that  d = 10  cm,  Ech  = 5.2  sv/cm,  and  r]  = 1 . 8x1  O'4  poise.  Then  Eq. 
(5.35)  may  be  written: 


V(volts)  = 


130500 


[2.088xl0'^B,n,  + D(|im)jt(sec) 


(5.36) 


Specializing  to  R=0.1//m  particles,  so  that  (Table  5.1)  B,n,  =6.949x10®,  we 


obtain 


V(volts)  = 


130500 

[1.451  + D(//m)]t(sec) 


(5.37) 
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Note  that  this  formula  also  can  apply  to  the  charging  section,  in  which  we 
assume  15600  volts  is  applied  across  10  cm,  and  the  particles  are  exposed 
for  one  second.  Provided  D is  less  than  6.91pm,  15600  volts  will  be  insufficient 
to  move  the  particles  across  the  channel  within  the  charging  time  of  one  second. 
Hence,  the  small-particle  charging  can  be  completed  before  the  particles  move 
to  the  walls  under  the  applied  charging  field. 


5.3.4  Maximum  Loading  17] 


max 


In  the  foregoing  equations,  the  mass  loading  17]  is  a "free"  parameter. 

It  cannot  be  arbitrarily  large,  however,  since  the  ions  would  eventually  become 
depleted  and  the  particles  would  compete  for  them.  It  is  of  interest  to  consider  the 
foregoing  equations  in  the  light  of  the  maximum  loading  17]  for  a given  ion 
supply  ( n+  ions/cm  ).  Under  conditions  where  all  available  ions  are  deposited 
on  particles,  and  assuming  that  most  of  the  charges  accumulate  on  the  large 
particles  by  conduction,  we  may  estimate  ]7]  by  equating  n = //,  n = 
Fe//onsat,  and  using  Eq.  (5.29)  for  7]q  and  Eq.  (5.30)  for  n . Thus,  solving 

for  V]=  77) 


max 


W] 


max 


= ^pe 

9 E 


-t-  = 1.933x1  (5.38) 

Ech  \F*  W \F2dV 


in  pm 


8 3 3 

where  we  have  used  n = 1.5x10  /cm  , p = 2 gm/cm  , and  E , = 5.2  sv/cm. 

t ch 

Using  this  value,  ^max*  hi  the  expressions  (5.32)  and  (5.34)  for 

and  T , respectively,  we  obtain  minimal  values  for  T and  T , and  a maximal 
R A R 

value  for  V from  (5.37).  The  minimal  values  of  T and  T_  are: 

A R 


/t.)  = (4^e2  n B,n,  ) * = 3.31  sec 

(for  R = 0.  lpm) 


min 


(5.39) 


(TR ) (=*■)  = 4* 808  )sec 

'Rtnin  en+EchW  \D(pm )/ 


(5.40) 


and  Vmax<27170  volts  (for  t=3. 31  sec  and  D = 0).  It  is  of  interest  to  note  that 
(*^)  -n  is  independent  of  the  distribution  of  Ec^. 
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5.3.5  Design  Criteria 


In  applying  the  above  formulas  to  the  design  of  an  agglomerator  it  seems 

reasonable  to  assume  maximum  loading  so  that  TA-  . This  minimizes 

max 

the  agglomeration  time  T . It  also  minimizes  the  self- repuls  ion  time  T 
(which  is  also  approximately  the  self-agglomeration  time  for  the  large  particles). 
The  following  relations  are  desirable,  namely,  that: 

(a)  rA  be  smaller  than  Tr;  then  agglomeration  can  occur  before  self- repulsion 
(and  before  self-agglomeration  of  the  large  particles). 

(b)  be  of  the  order  of  a second  (residence  time  in  the  agglomerator  section). 

(c)  V be  large  enough  to  bring  the  particle  clouds  together  and  mix  them  before 
they  disperse  by  self -repuls ion. 

Ech  *the  charSin£  ^eld)  be  small  enough  to  allow  small  particles  to  charge 
for  a second  or  so  before  they  are  pushed  to  the  walls  by  the  charging  field 
(charging  section). 

In  addition,  assuming  that  F2  /£>  is  not  too  small,  depending  on  the  dis- 
tribution, one  may  design  the  system  based  on  B,n,  for  particle  radius  anywhere 
in  the  approximate  range  R = 0.1/xm  and  R = 1.0pim,  since  B,n,  is  roughly  constant. 
Then,  since  B,n,  increases  rapidly  below  R - 0.1/xm  (due  mostly  to  the  Cunningham 
slip  factor),  the  agglomeration  times  7^  for  the  smaller  particles  become  smaller 

(inversely  proportional  to  Bjnj  ),  while  the  self -repuls ion  time  T remains  unaf- 

R 

fected.  In  other  words,  if  the  agglomeration  works  well  for  particles  of  radius, 
say,  R = 0.1pm,  it  will  work  even  better  for  smaller  particles.  As  Zebel  says 
(Ref.  5.6,  p 35):  ".  . . small  particles  are  very  rapidly  'eaten  up'  by  large  ones." 

5.4  Application  to  Typical  Size  Distributions 

We  may  apply  the  above  theory  to  some  sample  actual  size  distributions, 
in  order  to  Bee  how  these  distributions  affect  the  operation  of  the  system  (r  , 

V and  V).  We  have  chosen  the  first  7 distributions  from  Table  2.1  of  Sec.  2. 

(the  remaining  distributions  are  similar  to  one  of  the  7 chosen).  In  Table  5.2 
we  have  labelled  the  distributions  as  #1,  #2,  etc.,  labelling  the  columns  across 
the  top.  The  horizontal  entries  denote  the  following: 
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= mass  median  diameter 


= geometric  standard  deviation 
= geometric -mean  diameter  (count  median) 


= count -mean  diameter 


qz  = ratio  of  count-mean-cube  diameter  to  count-mean- squar  e diameter 


F2/D 


= fraction  of  large  particles 

= used  in  formula  for  (7V,) 

n mu 


= used  in  formulas  for  T and  17\ 


DW/F2 


= maximum  mass  loading 


7"  c)7]  - x agglomeration  time 

A 

T 7/i  = JfJ  x repul  sion  time 

K 

= minimum  agglomeration  time  ( = ^riax) 

(Tp)min  = minimum  repulsion  time  (ty  = ^max) 

IT  ) . I(tA  = ratio  of  minimum  times 

' R/min/  \ A/min 

V.  = V(volts)  to  mix  particles,  evaluated  with  t = (TV  ) 

A ' A/min 

V = V(volts)  to  mix  particles,  evaluated  with  t = (tr) 

K \ *Vmin 

We  see  that  IT.  . = 3.  31  3 sec  is  a r easonable  value  (for  the  time  of 
' A'min 

residence  in  the  agglomerator  section  such  that  significant  agglomeration  occurs) 
and  is  the  same  for  all  distributions.  The  voltages  V required  to  bring  the  par- 

x\ 

tides  together,  that  is,  to  mix  them  before  they  disperse  are,  except  for 
Case  #4,  under  about  23  kilovolts  (a  reasonable  maximum  voltage)  for  all  dis- 
tributions. For  Cases  #3,  #4,  #b,  and  #7,  (tr)  is  less  than  (t.)  . This 

means  that  once  the  particles  are  mixed,  there  will  be  some  agglomeration  of 
the  large  particles  amongst  themselves.  This  should  not  be  serious  for  Cases 
#3  and  #6. 


Table  5.2  Effects  of  Distribution 

(Distributions  from  Table  2.1) 


Dm/2("m| 


D (/Jm) 
g 


D (//m) 
\D3/D2)(, 


Fz 



ID3 

— = (//m) 
Ft  pi 

7n 

'max 

(g/cm3) 

ta^ 

(sec-p/cm3 

(sec-f^cm3) 

ftLn 

(sec) 

(sec) 

/r  ) 

V R mim, 


A 

(volts) 

VR 

( volt  8 ) 


1.8 

2.3 

4.0 

0.34 

2.7 

Dg  = 0.9 

2.2 

3.5 

1.7 

5.4 

2.4 

~3 

• 

0.2788 

0.02075 

1.7187 

6. 70 1x10 '5 

0.2709 

0.9 

0.3804 

0.04548 

1.9785 

2. 778xl0-4 

0.3974 

1.646 

1.319 

1.0495 

3.475 

0.0820 

1.8406 

18.38 

0.6633 

0.03525 

1.00 

F.  = 2.527 
1 xlO-5 

7.91xl0*6 

0.6356  ' 

0.9149 

1.7437 

0.7751 

0.5054 

0.02847 

1.5994 

0.5558 

1.989 

29.77 

3.475 

10367. 

2.896 

20.09 

3.845x10 

-5 

5.755x10 

6.717x10*  ^ 

0.02004 

-6 

5.598x10 

3.883x10 

1.274x10 

1.907xl0"4 

2.226x10*  5 

0. 0664 

1.855xl0"5 

1.287xl0-4 

3.223xl0'5 

2.145xl0-4 

1.633xlO-5 

2.744x10“ 

4.305xl0*5 

1.0379X10*4 

3.313 

3.313 

3.313 

3.313 

3.313 

3.313 

8.384 

3.727 

2.430 

0.1368 

7.690 

2.672 

2.530 

1.125 

0.7335 

0.04128 

2.  321 

0.8065 

21506 

26319 

11485 

27139 

21309 

12718 

8.341 

9.288 

14.28 

1.00 
• =4.41 
1 x!0~ 


3.313 

0.5172 


8498  23396  15659  657243  9180 
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Case  H 4 seems  to  be  highly  singular.  For  example: 

(1 ) is  less  than  one  Angstrom  unit 

(2)  The  fraction  of  large  particles  (F2)  is  excessively  small,  less  than  10-5. 
One  may  suspect  that  the  size-distribution  for  Case  *4  is  poorly  represent- 
ed by  a log-normal  distribution;  therefore,  the  small  ratio  (0.04)  of  r epulsion-to- 
aggl omeration  times,  and  the  excessive  voltage  (657  kilovolts)  to  bring  particles 
together,  may  be  unrealistic  for  this  case. 

In  summary,  of  the  7 cases  studied,  all  except  H and  #7  appear  to  be  rea- 
sonably amenable  to  agglomeration;  these  exceptional  cases  might  also  turn  out 
to  be  amenable  if  their  distributions  could  be  more  accurately  characterized. 

6.  EVALUATION 
6.1  Brief  Summary 

The  filtration  enhancement  system  considered  here  - which  consists  of 
charger  and  agglomerator  sections  (Figs.  2.1  and  2.2)  - and  the  scientific  areas 
requiring  investigation,  have  been  discussed  in  detail  in  Section  2.  The  basic 
objective  is  to  agglomerate  fine  particles  onto  larger  particles,  hence  improv- 
ing  prefilter  efficiency.  Sections  3-5  contain  the  initial  results  obtained  on 
each  of  the  areas  of  investigation,  which  can  be  summarized  briefly  as  follow  i. 

n r 

Results  on  the  ionization  distribution  due  to  the  charger  section  Kr 

beta  source  are  presented  in  Section  3.  The  beta  source  (Fig.  3.1)  was  fab- 
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ricated,  filled  with  0.1  atm  Kr  gas,  and  the  induced  ionization  around  the 
source  measured  in  a geometry  (Fig.  3.3)  similar  to  the  charger  section  (neg- 
lecting the  collimators).  Results  for  the  specific  ionization  are  given  in  Fig. 

3.4.  Also  given  in  that  figure  are  analytical  results  obtained  neglecting  all 
scattering  effects.  The  reasonably  good  agreement  suggests  that  for  calcula- 
tions where  only  the  energy  deposition  is  of  interest  (as  here),  rather  than  the 

electron  distribution  function,  this  is  a good  assumption.  Based  on  these  re- 
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suits  a source  strength  of  about  500  Ci  Kr  was  selected  as  a basis  for  further 
consideration,  leading  to  an  ion  pair  production  rate  near  10^/(cm^-sec)  at 


„ about  10  cm  radially  from  the  source. 
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The  ion  density  distribution  in  the  charger  was  investigated  in  Section  4, 
with  approximate  results  given  in  Figs.  4.2  and  4.3.  It  is  clear  from  these  fig- 
ures that  regions  of  well  - separated  positive  and  negative  ion  distributions  can 
be  maintained  near  the  oppositely  charged  plates  (as  in  Fig.  2.2)  with  applied 
potentials  of  <20  kV.  Thus,  the  air  to  be  filtered  is  separated  into  two  streams, 
one  passing  through  a region  of  positive  ions,  the  other  through  a region  of  neg- 
ative ions.  The  charging  process  for  the  particulates,  and  the  subsequent  ag- 
glomeration of  the  oppositely  charged  particulates  in  the  presence  of  the  reversed 
electric  fields  (Fig.  2.2),  are  considered  in  Section  5. 

The  operating  feasibility  of  induced  agglomeration  has  been  assessed  on 


the  basis  of  theoretically  derived  time-scales  for  charging,  agglomeration,  and 
space-charge  repulsion.  For  logical  choices  of  the  parameters,  we  obtain 

reasonable  agglomeration  times  (residence  time  in  the  agglomeration  section 
such  that  significant  agglomeration  occurs,  and  therefore  successful  operation) 
of  the  order  of  seconds.  In  the  agglomeration  process,  the  small  particles  are 
consumed  by  the  larger  particles.  With  plate  voltages  on  the  order  of  20  kV 
the  particles  can  be  oppositely  charged  in  the  two  halves  of  the  gas  stream,  and 
then  the  oppositely- char ged  clouds  of  particles  can  be  brought  together  (mixed) 
in  the  agglomerator  section  - in  a time  shorter  fhan  that  in  which  they  would 
"explode"  or  disper  se  under  their  space-char  ge  r epulsion  forces.  Thus,  dis- 
persion is  avoided,  as  well  as  the  self-agglomeration  of  the  large  charged  par- 
ticles with  one  another  (which  takes  place  in  the  same  order  of  time).  The 
dispersal  time,  and  therefore  the  voltage  required,  depends  on  the  particle 
size  distribution,  but  the  agglomeration  time  does  not.  We  have  considered 
0. 1 -jim-radius  particles,  and  have  argued  that  (a)  smaller  particles  than  these 
would  be  consumed  even  faster  by  the  large  particles,  and  that  (b)  particles  with 
radii  in  the  range  0.  l/<m  to  l.Ofjm  would  be  agglomerated  onto  the  "supermicron" 
(above  Igm)  particles  about  equally  rapidly.  Thus,  it  is  expected  that  the  size 
distribution  would  shift  significantly  toward  the  supermicron  sizes,  which  is 
the  desired  goal. 
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Asa  consequence,  we  believe  that  the  radiation  induced  particulate 
agglomeration  system  considered  here  could,  in  fact,  lead  to  significant  en- 
hancement of  the  prefilter  efficiency,  and  hence  overall  improvements  in 
economy  in  operation  of  nuclear  facility  filtration  systems. 

6.2  Comparison. of  Results  with  Some  Related  Work 

Heinsohn  et  al.  (Ref.  2.1)  have  carried  out  an  investigation  of  possible 
, 60  . . 

use  of  Co  gamma  radiation  for  production  of  ionization  in  a particulate  charg- 
ing system.  Their  objective  was  to  remove  the  particles,  after  being  charged, 
by  electrostatic  precipitation  - quite  a different  principle  from  the  agglomera- 
tion technique  investigated  here.  Yet,  some  of  their  results  on  ionization  are 

quite  useful  for  comparison.  Their  experimental  system  used  two  75  Ci  sources 
60 

of  Co;  the  ionization  was  produced  by  Compton  and  photoelectric  electrons 
ejected  in  their  experimental  housing  by  the  ~1  MeV  ^°Co  gamma  rays.  Their 
maximum  measured  specific  ionization  was  about 

9.3  73 

4.8x10  [ions/(cm  -sec)]/(2x75  Ci)«3xl0  ions/(cm  -sec-Ciy's)  (6.1) 


Note  that  this  is  several  times  less  than  the  specific  activity  in  the  10-20  cm  re- 
gion of  Fig.  3.4.  Furthermore,  the  radiation  hazard  external  to  the  apparatus 

85 

that  is  associated  with  Kr,  when  measured  per  curie,  is  far  less  than  that  due 
60  85 

to  Co.  The  Kr  gamma  ray  (~0.5MeV)  is  emitted  in  only  0.4%  of  the  decays, 

and  the  efficiency  of  conversion  of  the  99.6%  emitted  beta  energy  into  bremsstrah- 

lung  (£670keV,  avg~3 00  keV)  is  less  than  1%.  Thus,  on  a radiation  hazard  basis 

Kr  is  less  than  about  1%  of  a similar  source  strength  ^Co  source.  We  concur 

with  Heinsohn  et  al.  (Ref.  2.1),  who  note  in  their  conclusion  that  since  it  is,  in 

fact,  electrons  which  produce  the  ionization,  they  can  more  easily  be  provided 

85 

by  beta  than  gamma  radiation.  We  note,  for  example,  that  the  500  Ci  Kr  beta 

11  3 

source  found  here  to  provide  the  ~10  ions/(cm  -sec)  needed  ionization  rate, 
is  equivalent  to  only  5 Ci  of  external  (to  the  apparatus)  electromagnetic  radiation 
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( y's  and  bremsstrahlung).  It  is  also  of  interest  to  note  that  the  resulting  ion 
densities  are  similar  to  those  typically  found  in  electrostatic  precipitators 
(Refs.  5.2  and  6.1),  but  that  our  method  of  production  of  the  ionization  requires 
much  less  power. 

Regarding  electrically-induced  agglomeration  studies  by  Melcher  and 
Sachar  (Ref.  6.  £),  they  did  not  consider  radii  below  one  micrometer,  and  hence 
did  not  consider  fine  particles.  They  artificially  chose  two  groups,  a 1.0-  |im- 
radius  group  to  represent  all  the  "small"  particles,  and  a 1 0. 0- pm- radius  group 
to  represent  all  the  "large"  particles  in  a distribution;  whereas  we  consider  a 
series  of  entire  log-normal  distributions,  all  based  on  actual  nuclear  facility 
particle  distributions.  (Our  small-particle  fraction  as  a function  of  particle 
radius  is  assumed  to  be  given  by  the  integral  over  *he  log-normal  distribution, 
from  zero  radius  to  the  radius  of  interest.  This  is  based  on  the  analytical 
fall-off  behavior  of  the  log-normal  distribution  at  small  radii,  and  may  be  con- 
sidered a natural  way  of  characterizing  the  small-particle  fraction.)  In  Ref. 

6.2  it  was  assumed  that  the"small"and"lar ge"oppositely-char ged  particles  were 
somehow  instantaneously  charged  and  mixed.  This  seems  somewhat  artificial 
as  opposed  to  our  method  of  charging  and  of  bringing  the  particles  together, 
under  controlled  conditions.  Also,  in  Ref.  6.2  a special  ratio  of  dust  loading  dis- 
tributed between  the  'bmall"and"lar  ge"particles  is  postulated  in  order  to  achieve 
efficient  operation;  it  is  difficult  to  see  how  this  ratio  can  be  controlled  in  prac- 
tice. In  our  method,  the  total  dust  loading  is  the  (more  practical)  parameter 
affecting  the  rate  of  agglomeration.  The  values  of  dust  loading  assumed  in 
Table  5.2  leading  to  successful  operation  are  typically  found  in  industrial  elec- 
trostatic precipitators  (Ref.  6.1). 
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7.  RECOMMENDATIONS 

Based  on  the  above  evaluation,  we  believe  that  experimental  and  sup- 
porting analytical  work  should  be  carried  out  to  make  preliminary  tests  of  the 
charger  section,  as  the  first  part  of  Phase  II.  In  particular,  the  following 
tasks  can  be  defined: 

1)  Design  experimental  model  charger  section  and  fabricate. 

2)  Measure  ion  density  distribution  in  charger  section  (Langmuir  probe 
or  similar  instrument). 

3)  Modify  theory  based  on  Phase  I results  and  apply  to  obtain  better 
estimates  of  particulate  charging  and  agglomeration. 

4)  Measure  actual  state  of  charging  on  simple  (0.3/xm  monodisper se) 
particulates  (mobility  analyzer  or  similar  instrument). 

5)  Analyze  results,  write  report  and  suggest  agglomeration  studies 
if  charging  results  justify  it. 

It  is  expected  that  item  4)  would  be  carried  out  by  use  of  laboratory 
facilities  at  the  Harvard  School  of  Public  Health,  under  the  overall  supervision 
of  Dr.  Melvin  W.  First.  It  is  intended  that  the  minimum  source  intensity  will 
be  used  that  is  sufficient  to  demonstrate  the  charging  principles.  Maximum  use 
will  be  made,  in  the  experiments,  of  lower  than  actual  flow  rates,  since  the  flow 
rate  is  a parameter  that  occurs  in  the  theory.  Once  the  system  concept  has  been 
verified  for  low  source  strengths  and  flow  rates,  results  can  be  predicted  with 
confidence  for  realistic  system  parameters. 


53 


REFERENCES 


2.2 


3.1 

3.2 

3.3 

5.1 

5.2 

5.3 


5.4 

5.5 


t 


6.1 


6.2 


Heinsohn,  R.  J.,  Levine,  S.  H.,  Fjeld,  R.  A.  and  Malamud,  G.  W., 
"Radiation  Charging:  A Novel  Way  to  Electrically  Charge  Fine  Particles," 
409-46  in  Symposium  on  Electrostatic  Precipitators  for  the  Control  of 
Fine  Particles,  EPA-650/2-75-016,  C.  E.  Feazel,  Editor  (1975). 

Nelson,  G.^O.,  Miller,  B.  H. , Richards,  C.  P.,  and  Bierman,  A.  H., 
Enhanced  Filtration,  July-Sept.  1976  Progress  Report,  Lawrence  Liver- 
more Laboratory  Report  No.  UCID-16949-76-3  (March,  1977). 

Varma,  M.  N.  and  Baum,  J.  W.,  New  "wall-less"  ionization  chamber  for 
air  dose  measurement,  Rev.  Sci.  Instru.  47  112-114  (1976). 

Hogan,  O.  H.,  Zigman,  P.  E.,  and  Machin,  J.  L.,  Beta  Spectra  II.  Spectra 
of  Individual  Negatron  Emitters,  USNRDL-TR-802,  (Dec.  1964).  AD  455961. 

Berger,  M.  J.  and  Seltzer,  S.  M.,  Tables  of  Energy  Losses  and  Ranges  of 
Electrons  and  Positrons,  NASA  SP-3012  (1964). 

Dennis,  R.,  editor.  Handbook  on  Aerosols,  ERDA  Report  TID-26608  (1976). 

White,  H.  J.,  Industrial  Electrostatic  Precipitation,  Addison-Wesley, 
Reading,  Mass.  (1963). 

Parker,  L.  W.,  "Field  Charging  Theory  for  Electrostatic  Precipitation  of 
Ultrafine  Particles,"  68th  APCA  Annual  Meeting,  Boston  (1975);  also, 

Liu,  B.  Y.  H.  and  Yeh,  H-C,  J.  Appl.  Phys.  39  1 396  (1968). 

Davies,  C.  N. , Air  Filtration,  Academic  Press,  New  York  (1973). 

Parker,  L.  W.,  in  Precipitation  Scavenging  (1 974),  R.  G.  Semonin  and 
R.  W.  Beadle,  editors,  ERDA  Symposium  Series  No.  41  (CONF- 741 003), 
pp.  276-298  (1977). 

Zebel,  G.,  in  Aerosol  Science,  C.  N.  Davies,  editor,  Academic  Press, 

New  York,  pp.  31-58  (1966). 

Oglesby,  R.  and  Nichols,  G.  B.,  Southern  Research  Institute  Report, 

A Manual  of  Electrostatic  Precipitator  Technology,  Part  I,  Fundamentals, 
PB-196380  (1970). 

Melcher,  J.  R.  and  Sachar,  K.  S. , Electrical  Induction  of  Particulate 
Agglomeration,  MIT  Report  PB-205188  (1971). 


I 


* 


54 


